function [cineq, ceq] = nonlcon(c)

    global delta sigma_l sigma_u SampleSize
% Nonlinear inequality constraints
% Note: we use a closed form equivalent of equations (4) and (5) on page
% 1314.
    cineq = [    
    0.95-normcdf(c(1),0,1)+normcdf(-c(2)-(sqrt(SampleSize)*(delta)/sigma_u),0,1);
    0.95-normcdf(c(2),0,1)+normcdf(-c(1)-(sqrt(SampleSize)*(delta)/sigma_l),0,1)];
 
% Nonlinear equality constraints.
    ceq = [];

end